Comparing pertinent effects of antiferromagnetic fluctuations in the two and three 
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We use the dynamical vertex approximation (DFA) with a Moriyaesque A correction for studying 
the impact of antiferromagnetic fluctuations on the spectral function of the Hubbard model in 
two and three dimensions. Our results show the suppression of the quasiparticle weight in three 
dimensions and dramatically stronger impact of spin fluctuations in two dimensions where the 
pseudogap is formed at low enough temperatures. Even in the presence of the Hubbard subbands, 
the origin of the pseudogap at weak-to-intermediate coupling is in the splitting of the quasiparticle 
peak. At stronger coupling (closer to the insulating phase) the splitting of Hubbard subbands is 
expected instead. The k-dependence of the self energy appears to be also much more pronounced 
in two dimensions as can be observed in the k-resolved DFA spectra, experimentally accessible by 
angular resolved photoemission spectroscopy in layered correlated systems. 

PACS numbers: 71.27.+a, 71.10.Fd 
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I. INTRODUCTION 

Since its formulation, -'^ the Hubbard model served as 
a minimal model for electronic correlations. Due to the 
complexity of electronic correlations, solving this model 
is however only possible in dimension d — 1 (exactly via 
the Bethe Ansatz^) and in the limit d — oc^''*''^ (where 
the mapping^ onto an Anderson impurity model allows 
for an accurate numerical solution^i^) . Of physical in- 
terest are however strongly correlated systems in d = 3, 
for modeling the Mott-Hubbard transition^ and (anti- 
)ferromagnetisniiii2ii^, and \n. d — 2 for describing the 
cupratesiS, where the role of the antiferromagnetic fluc- 
tuations in developing pseudogap structures and super- 
conductivity are at the center of attention. 

The aim of this paper is to study the difference be- 
tween the effect of antiferromagnetic fluctuations on the 
electronic properties in d = 2 and d — i. For weak cou- 
pling (small Coulomb interaction U), the perturbation 
theory, and its extensions, e.g. the fluctuation-exchange 
approximation (FLEX)^^, the two-particle self-consistent 
approximation (TPSC)^^, and the functional renormal- 
ization groupie are suitable methods for this purpose. 
In d = 3 antiferromagnetic fluctuations produce only 
quantitative changes of electronic spectrum, although the 
particle-hole excitations enhance the quasiparticle scat- 
tering rate when the temperature T is approaching the 
Neel temperature. In d = 2 there are divergences in the 
self-energy diagrams and the abovementioned approxi- 
mations predict pseudogap structures in the self-energy 
in the weak-coupling regime^^^ii^. These techniques are 
however not applicable at stronger coupling, since they 
do not describe strong quasiparticle renormalization due 
to the Mott physics. 

Since we are interested in intermediate-to-strong elec- 
tronic correlations, we need to take a different approach. 
Starting point is the by-now widely employed dynam- 



ical mean-field theory (DMFT).3'4-5 This method be- 
comes exact^ for d — > 00, and yields a major part of 
the electronic correlations, i.e., the local correlations. 
However, any non-local correlations are neglected and 
hence DMFT does not differentiate between the Hubbard 
model in two- and three dimensions. More precisely, only 
differences stemming from different shapes of the density 
of states (DOS) are taken into account, not those result- 
ing, e.g., from antiferromagnetic correlations since these 
correlations are by nature non-local. 

Hitherto, the focus of DMFT extensions has been on 
short-range correlations within a (finite) cluster instead 
of the single DMFT impurity site. These cluster exten- 
sions of DMFT^^ have been used for describing pseudo- 
gaps and superconductivity in the two-dimensional Hub- 
bard model. Due to numerical limitations, the inclusion 
of important long-range correlations and the application 
of this method in three dimensions or realistic multi- 
orbital calculations is however not possible, except for 
very small clusters with 0(2 4) sites. Also the 1/d 
expansion of DMFTi^ is restricted to short -range corre- 
lations, as is a recent perturbative extensioni^^ 

Hence, for including long-range correlations, the focus 
of the methodological development has shifted recently 
to diagrammatic extensions of DMFT such as the dy- 
namical vertex approximation (DrA)^i2ii2^iS2i and the 
dual fermion approach by Rubtsov et al^ Even before, 
Kuchinskii et al^ combined the local DMFT self en- 
ergy with the non-local contributions to self energy of 
the spin-fermion model, and included long-range correla- 
tions this way. Their procedure, however, does not rely 
on a rigorous diagrammatic derivation. 

To include long-range fluctuations in a diagrammatic 
way DFA considers the local vertex instead of the bare 
interaction. It includes DMFT but also long-range corre- 
lations beyond. Our understanding of the physics asso- 
ciated with such long-range correlation is typically based 
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on ladder diagrams, which are considered, e.g. by the 
abovementioned TPSC and FLEX approximations. For 
example, the ladder diagrams in the particle-hole channel 
yield antiferromagnetic fluctuations in the paramagnetic 
phase (paramagnons) and (anti-)ferromagnons in the or- 
dered state. It is natural to suppose that the contribu- 
tion of the corresponding fluctuations in the intermediate 
coupling regime can be described by the same kind of di- 
agrams albeit with the renormalized vertices. In DFA 
the local (frequency dependent) vertex is considered in- 
stead of the bare interaction. Therefore, this method 
reproduces the results of the weak-coupling approaches 
at small U but can treat spatial correlations also at inter- 
mediate coupling. Hence, DFA is well suited for studying 
antiferromagnetic fluctuations in strongly correlated sys- 
tems both for d — 2 and d = 3. 

The paper is organized as follows: In Section [TT] we 
reiterate the DFA approach in a formulation with the 
three-point (instead of the four-point) vertex functions 
which allows for a connection to the spin fermion model 
in Section Hm and for the analytical considerations on the 
DFA self energy in Section IIVI In Section |Vl we intro- 
duce a Moriyaesque A correction to the susceptibility to 
describe correctly the two-dimensional case. Results for 
three dimensions are presented in Section IVII and com- 
pared to those in two dimensions in Section [VIII Special 
emphasis to angular resolved spectra is given in Section 
IVIIII before we give a brief summary in Section IIXI 




FIG. 1: (Color online) Graphical representation of the con- 
tribution of bare Coulomb interaction (a) and spin (charge) 
fluctuations (b) to the self-energy in the DFA approach, Eq. 
([8]). Solid lines correspond to the electronic Green function 
Gk,i^, dashed line to the bare Hubbard interaction U , wiggly 
lines - to the spin (charge) susceptibility Xq^i the triangle 
corresponds to the interaction vertex 7^^'^ ^. 



can be expressed in terms of the particle-hole bubble 

Sioc(i^)]~^ is the Green function, and Yj\oc{v) the lo- 
cal self-energy. The spin (charge) irreducible local ver- 
tices F^j'^^^jj. are determined from the corresponding local 
problemSS. 

The result ^ accounts for the contribution of ladder 
diagrams to the self-energy in the two particle-hole chan- 
nels. Following Edwards and Herta^ it is convenient to 
pick out parts of these ladders, which are separated by 
the bare on-site Coulomb interaction U. This is achieved 
by considering the quantities 



s(c),q 



(4) 



W,q 



II. DYNAMICAL VERTEX APPROXIMATION 

Starting point of our considerations is the Hubbard 
model on a square or cubic lattice 



H 



^ .t 



(1) 



where t denotes the hopping amplitude between nearest- 
neighbors, U the Coulomb interaction, cl^{cia) creates 
(annihilates) an electron with spin a on site i; riia- = 
cl^Cia- In the following, we restrict ourselves to the para- 
magnetic phase with n — 1 electrons / site at a finite tem- 
perature T. 

The DFA result for the self-energy of the model ([T|) 
was derived in Ref. [13, see Eq. (16). For the purpose 
of the present paper this result for the self-energy can be 
written in the form 



such that x^'c'^q = {['^Xi^r^ T U}-^ with the upper 
(lower) sign for the spin (charge) susceptibility. The non- 
local spin (charge) susceptibility is then given by 



(5) 



and therefore 



provided to be a particle-hole irre- 



ducible susceptibility in the spin (charge) channel. Intro- 
ducing, similar to Ref. 28, the corresponding three-point 
vertex 7g(^) q of electron interaction with charge (spin) 
fluctuations, 



7s(c),q 



*(c),q' 



(6) 



the irreducible susceptibility (/)q,^j'' can be represented as 



7s(c),qX0q 



(7) 



^Xs,ci s.ii' 



Ac,q c,ir 



XoqLjv c,loc s.loc / 



In these notations, the result ^ can then be rewritten 
identically as 



(2) 



where the non-local spin (s) and charge (c) suscepti- 
bilities 



As(c), q ~ LlXOqo)/' "yu s(c), irJ 



(3) 



+ 3C/7r,qX^. + f/7c.qXq. 
^ / ^ XOqcj \ c,loc s,loc / 



c,q 



G 



k+q,z-'+LJ 



(8) 
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The first three terms in the square brackets correspond to 
the interaction of electrons via Hubbard on-site Coulomb 
interaction (without forming ph-bubbles, Fig. [T^), the 
next two terms correspond to electron interactions via 
charge- and spin- fluctuations (Fig. [T|d), the last term 
subtracts double counted local contribution. 



III. RELATION TO SPIN-FERMION MODELS 
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The contributions of bare Coulomb interaction and 
charge (spin) fluctuations to the self-energy ([8]) can be 
also obtained from the fcrmion-boson model with gener- 
ating functional 

Z = J D[cl^,Cka]DSq^^Dpq^^exp{-£[S, p,c]} 
C[S, p,c]= ^ (w„ - ev,)c\^Ck„ (9) 

~\~U ^ ^^ {Pqu)P—q,—uJ "I" ^quj^ — q,—uj) 
q,a; 

+ U ^ (7s,q)^''^4,l/,a'^'7'7'Ck+q,l/-fc^,CT'Sq,a; 

k,q,z-',cj,cj,(7' 

+iU ^ (7c,q)^^^4.l.,aCk+q,l.+^,<7Pq,.. 
k.q.l/.oj 

where „ is determined in the present approach ac- 
cording to the Eq. ([6]) and (Jaai are the Pauli matrices. 
The model © is similar to that derived from the Hub- 
bard model via Hubbard-Stratonovich transforniation3^, 
but it is explicitly spin symmetric and contains the non- 
local frequency dependent vertices lc(s) q' which account 
for the local- and short range-nonlocal fluctuations. 

Contrary to the earlier paramagnon theories^ and the 
spin-fermion mode l^^'^^ , where 7^ ^ = 1 and charge fluc- 
tuations are omitted (7^ q = 0), we have 7^^"^ q 7^ ^^'^ 
^ 1. The frequency dependence of the vertices l^sfc) q 
calculated in the present approach for two dimensions 
with Q = (tt, tt) is shown in Fig. [2] (in the three di- 
mensional case we observe qualitatively similar behav- 
ior). One can see, that both charge- and spin vertices 
have a strong frequency dependence and approach unity 
only in the high-frequency limit. While in the weak- 
coupling regime U = D = At both vertices are sup- 
pressed at small frequencies [which is the consequence 
of the particle-particle (Kanamori) screening], closer to 
the DMFT Mott transition (at U = 2D = 8t) the spin 
vertex at small frequencies is enhanced. This behavior is 
similar to that observed in Ref.ll^for the three-frequency 
(four-point) vertex in the three dimensional case. 

Hence, the spin-fermion theory, which was heuristically 
added to the DMFT self-energy before, is included in a 
more systematic and consistent way in DFA, which also 
accounts for the corrections to the electron-paramagnon 
vertex. The susceptibility Xq uj which is determined phe- 
nomenologically in the spin-fermion model is obtained 



FIG. 2: (Color online) Frequency dependence of the spin and 
charge three-point vertex Eq. ([SJ at (7 = 1, /3 = l/T = 15 
(left) and J7 = 2, /3 = 10 (right), = 0, at the antiferromag- 
netic wave vector Q = (tt, tt). All energies are in units of half 
the effective bandwidth D = it. 



in our approach by dressing the bare propagator 1/U of 
charge- and spin fields by particle- hole bubbles, which re- 
produces the results ([5]) and (O of the previous Section. 

Using the model (O one can also calculate the leading 
order non-local correction to the three-point vertices due 
to fermion-boson interaction, 



I s,k,q 



' S.Q ' 



Itu y 

2 ^ 

1^1, qi 



/ s,q 



[2- 



/s.qi 



/c,qi 



'U-y"'^^ y" -I- C/7''"! 



G 



k-(-qj,l^-|-aJi 



I c,k,q 



7r,q + 



Itu y 
2 ^ 



loc, 



(10) 



/s,q L ' 



'c,qi 



+3C/7r,q^ 



loc. 



hqi,y+wi 

hqi+q.l'-ftxJl-l-ixJ " ^"-^^J (-'--'-) 

where loc stands for the subtraction of the local terms 
already included in 7^ q. The non-local corrections to 
the self-energy and vertex can be then treated self- 
consistently by substituting them into Eq. ([7]). This pro- 
vides an alternative simpler way of self-consistent treat- 
ment instead of the more complicated parquet approach 
discussed in Ref. An even simpler way to go beyond 
a non-self consistent treatment of the DFA equations is 
considered in Sect. V. 



IV. ANALYTIC APPROXIMATION FOR THE 
DFA SELF ENERGY 

Similarly to the weak-coupling approach^^, in the two 
dimensional case the self-energy can be obtained approx- 
imately analytically. In this case the susceptibility x^i^ 
is strongly enhanced at = and q « Q = (tt, tt), and 
can be represented in the form 



XqO 



A 



(12) 



(q-Q)2+r' 

where = ^/(l " t/0Qo) with A - (V^q^'^oU^Q being 
the (squared) inverse spin fluctuation correlation length. 
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Since the corresponding momentum sum in the Eq. ([8]) 
over q is logarithmically divergent at ^ — > oo, we can 
approximately retain ourselves to only the zero bosonic 
Matsubara frequency term in the spin-fluctuation contri- 
bution and put q w Q in all the factors except Xqo to 
obtain 



(13) 



where = (3^2/2) ^ Xq.o- 



To study the frequency dependence of the self-energy 

s,Q 



qualitatively, we first consider 7^ q = 1 and choose 



the local self-energy in the form (see, e.g. Ref. |3S 

S,o,(^) = (1 - ^)(AL/4)/(^ - AL«:/(4z.)) (14) 

where Aioc ~ [/ is the size of the Hubbard gap and n 
measures the relative weight of the quasiparticle peak 
(QP) with respect to the Hubbard subbands (k = at 
the Mott transition and k = 1 for [/ -> 0). The Eq. 
allows to reproduce the three-peak structure of the self- 
energy, observed in the numerical solution of the single- 
impurity Anderson model, supplemented by the DMFT 
self-consistent condition. 

The evolution of the spectral properties calculated 
with the self-energies and (I14p with changing k for 
Aloe = 1 and A = 0.1 is shown in Fig. [3] (we suppose 
that the vector k is located at the Fermi surface and 
Ek+Q = due to nesting). One can see that at small k, 
i.e. in the vicinity of the Mott transition one finds split- 
ting of Hubbard subbands, while the QP remains unsplit 
(Fig. 3a, b). In the narrow region of larger k the QP 
is split in two peaks, and the splitting of the Hubbard 
subbands remain visible (Fig. 3c). At intermediate val- 
ues of K we find only splitting of the QP peak, the two 
other peaks corresponding to the Hubbard subbands are 
present (Fig. 3d,3e). Finally, in the weak coupling limit 
K = 1 we reproduce the two-peak pseudogap, discussed in 
Refs. [l^fisl (Fig. 3f). In a more general case of 7^q 7^ 1 

we expect a pseudogap of the size ~ Al^yf^^y^^ in the 
weak coupling regime at small enough temperatures and 
more complicated structures at strong {7; see our numer- 
ical results below. 



MORIYAESQUE A CORRECTION FOR THE 
VERTEX 



The local approximation for the particle-hole irre- 
ducible vertex, considered in Section II, is however not 
exact. In particular, the magnetic transition tempera- 
ture remains equal to its value in DMFT, and therefore 
it is overestimated in both three- and two dimensions. In 
the latter case T/v would remain finite, contrary to the 
Mermin- Wagner theorem. 

In the DFA framework a reduction of Tjy would nat- 
urally arise from a self-consistent solution of the DFA 
equations. An alternative (simpler) way to fulfill the 





FIG. 3: (Color online) The spectral functions in d = 2 as 
obtained from the approximate self-energies including local 
(dashed lines, Eq. (|14|l ') and non-local (solid lines, Eq. (I13|) ) 
fluctuations for «: = (a), 0.1 (b), 0.3 (c), 0.5 (d), 0.9 (e), and 
1.0 (f). 



Mermin- Wagner theorem in 2D (and to reduce the tran- 
sition temperature in three dimensions) is to introduce a 
correction to the susceptibility similar to the Moriya the- 
ory of weak itinerant magnets^S. To this end, we replace 



(15) 



Formally the r.h.s. of Eq. is exact for some (un- 

known) Aqi^; in the following we assume Aq^^ ~ Aqo = A 
since static fluctuations with momentum Q predominate 
near the magnetic instability. Instead of determining (as 
it was done in Moriya theory) A from the fluctuation cor- 
rection to the free energy, which is rather cumbersome in 
the present approach, we (similar to TPSC) impose the 
fulfillment of the sumrule 



— ImSk.i. = U'^n{l - n/2)/2. 

TT 



This also implies 



ReSk,!/ 



f/2n(l - n/2) 
2^^ 



(16) 



(17) 



for v ^ D, according to the Kramers-Kronig relation. 
The latter asymptotic behavior may be very important 
to obtain the correct Fermi surface in the non-half-filled 
case, but should be fulfilled also in the half-filled case to 
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FIG. 4: (Color online) DFA self-energy on the Matsubara 
axis calculated with and without Moriya A correction for two 
different points of the Fermi surface in the two dimensional 
Hubbard model {at U = D = 4t, [3 = 1/T ^ 17); also shown 
is the DMFT self energy for comparison. Notice that, with- 
out introducing the Moriya A correction, one always observes 
a deviation of the high-frequency E(k, i!/„) from the correct 
asymptotic behavior ~ C/^n(l — ^)/(2ii/„) = U'^/{Aiu„) which 
is consistent with the self energy sum rule (see text). 



obtain correct spectral functions. It is obviously violated 
in standard spin-fermion (also paramagnon) approaches 
in two dimensions, where the Neel temperature (T^r) is 
finite without the A correction and the l.h.s. of Eqs. ITBl) 
and pT)) are divergent at T — > Tm- 

The frequency dependence of the self-energy at the 
imaginary axis for the two-dimensional Hubbard model 
{U = D — At), calculated with and without A correc- 
tion is compared in Fig. 2) The A correction removes 
the divergence of the l.h.s. of Eqs. (|16p and (fT7|) at 
T T^^^^ and leads to the correct asymptotic behav- 
ior at large Vn. Without A-correction (or, alternatively, a 
self-consistent solution of the DFA equations) spin fluc- 
tuations and their pertinent effect on the self energy are 
overestimated. This is because the spin fluctuations re- 
sult in a reduced metalicity which in a second DFA iter- 
ation, i.e., the recalculation of the local vertex with the 
less metallic Green function as an input^^, would reduce 
the spin fluctuations. 

In two dimensions the sumrules (fTH)) and (flT)) can be 
fulfilled at all positive temperatures, and the actual tran- 
sition temperature is zero, as required by the Mermin- 
Wagner theorem. As one can see from Eqs. ([T^ and 
(fT3)) . the correlation length ^ in two dimensions is expo- 
nentially divergent (with A-correction): 



e(xexp(6/r), 

where the coefficient b in the exponent is proportional 
to U . This is evidently confirmed also by our numerical 
results shown in Fig. [5l where we have reported the val- 
ues of the inverse of the spin susceptibility at Q = (tt, tt) 
calculated with the inclusion of the A-correction: The 
exponential divergence of ^ for T ^ is directly re- 



FIG. 5: Temperature dependence of the A-corrected inverse 
antiferromagnetic susceptibility in two dimensions (triangles) 
for U = D = it. The data display an exponential temperature 
dependence, consistent with the expected behavior of ^ (see 
text). The DMFT Neel temperature corresponding to this set 
of parameters is marked with an arrow. 



fleeted in an analogous behavior of the spin susceptibility 
(Xq.o - M^, see Eq. ([12])) at T ^ 0. 

In three dimensions, on the other hand, the sumrules 
(HH) and dni) with (xq^o)"^ + a > can be fulfilled only 
down to a certain temperature T^^^, which is reduced in 
comparison with T^^^^ and determines the phase tran- 
sition temperature in the DFA approach. 



VI. RESULTS FOR THE HUBBARD MODEL IN 
THREE DIMENSIONS 

Let us turn to the results for the self-energy and spec- 
tral functions which are obtained applying the Moriya 
A correction to the vertex of the DFA for the three di- 
mensional system (the analytical continuation to the real 
axis ivn — > w was done using the Fade algorithm). In this 
case, as mentioned above, the A correction is expected to 
result in small -and only quantitative- changes of the final 
DFA results, because in d = 3 (where the antiferromag- 
netic long-range order survives at finite temperatures) 
the A correction produces just a moderate reduction of 
the Neel temperature w.r.t. the DMFT value. 

Our results, shown in Fig. [HI clearly confirm this ex- 
pectation. Specifically, we analyze the case, already con- 
sidered in our previous study Ref. [13, i.e., the three di- 
mensional Hubbard model with U = 1.5 (in the units 
of half the variance of the non-interacting DOS, being 
D = 2V6t for d = 3), and /3 ^ 11.2 (in units of 
l/D), which corresponds to a temperature slightly above 
the DMFT Neel temperature (T^^^^), but appreciably 
higher than the three-dimensional T^^^ with A correc- 
tion (an estimate of the A— reduced Neel temperature 
gives /3°r^ = ~ 16.5). In this situation, as 

noticed in Ref. and shown in Fig. [S] (first row), the 
standard DFA results display a sizable renormalization of 
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FIG. 6: (Color online) DMFT self-energies and spectral func- 
tions (grey dashed line) at Icf = (7r/2, 7r/2, 7r/2) for the Hub- 
bard model in d = 3 at [/ = 1.5D [D = 2y/Qt) and (3 = 11.2 
(i.e., slightly above T^^^^) are compared with the corre- 
sponding DFA results with (lower row; solid blue line) and 
without (upper row; black dotted line) Moriya A correction. 
Note that (i) the non-local fluctuations modify only quanti- 
tatively the shape of the QP, but no pseudogap appears, and 
(ii) non-local correlation effects are further reduced by the 
inclusion of the Moriya A correction. 



the quasiparticle (QP) peak present in the DMFT spec- 
trum. However, no qualitative change in the nature of 
the spectral functions can be observed. The inclusion of 
the Moriya A correction, as shown in the second row of 
Fig. [6l reduces the renormalization effects due to non- 
local correlations: both the real and the imaginary part 
of the DFA self-energy at low frequency get very close 
to the DMFT values, and, obviously, the same happens 
to the QP peak in A(k,w). This result is easily under- 
stood in terms of the reduction of Tjv determined by the 
Moriya corrections, since the enhanced distance to the 
second-order antiferromagnetic transition at T/v leads to 
a reduction of the spin-fluctuation and corrections to the 
DMFT self-energy. If we reduce the temperature towards 
the DFA Neel temperature, antiferromagnetic spin fluc- 
tuations become strong again, and as shown in Fig. [3 
we indeed find results which are qualitatively similar to 
those without A correction (first row of FiglH]). In partic- 
ular, in both figures the quasiparticle weight is smaller 
in DFA than in DMFT in agreement with the expected 
effect of antiferromagnetic fluctuations. 

Summing up the results for the isotropic three dimen- 
sional system, we emphasize that the principal conse- 
quence of the inclusion of the Moriya A correction is a 
shift of the region with appreciable non-local correlation 
effects (i.e., the region where the DFA spectra substan- 
tially differ from DMFT) to lower temperatures, i.e., to 
the proximity of the "new" line of the antiferromagnetic 
phase transition. Our result demonstrates that for d = 3 
-with or without lambda correction- the extension of the 
region characterized by relevant non-local correlations is 
relatively small even for intermediate values of the inter- 
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FIG. 7: (Color online) DMFT self-energies and spectral func- 
tions at kF = (7r/2, 7r/2, 7r/2) for the Hubbard model in d = 3 
at 17 = 1.5D {D = 2y/&t) and /? = 15 compared with the cor- 
responding DFA ones with A correction. Lowering T = l/f3 
towards T^^"^ , qualitatively similar results as without A cor- 
rection at higher T (upper row of Fig. [6] /3 = 11.2) are ob- 
tained. 



actions. This indicates, hence, that for d = 3 DMFT 
represents indeed a good approximation, except for the 
region close to the antiferromagnetic transition. 



VII. RESULTS FOR THE HUBBARD MODEL 
IN TWO DIMENSIONS 

The effects of non-local correlations are -as one can 
imagine- much more dramatic for a two-dimensional sys- 
tem. It is easy to figure out that the divergence of the lad- 
der diagrams in the spin channel leads to huge non-local 
corrections in the DFA self-energy, which can differ also 
qualitatively from the DMFT one. At the same time, one 
should expect that in two dimensions the non-local corre- 
lation effects could be sensibly overestimated by the DFA 
without the inclusion of the Moriya A correction. As we 
have discussed in Section |Vl these corrections are essen- 
tial to fulfill the Mermin- Wagner theorem, pushing the 
Neel temperature from the DMFT value down to zero. 
Hence, for any finite temperature the antiferromagnetic 
fluctuations are reduced. The effects of the divergence of 
the spin ladder diagrams are also to some extent atten- 
uated in the formula for the DFA self-energy, because of 
the extra dimension gained at T = due to the trans- 
formation of the Matsubara summation to a frequency 
integral on the r.h.s. of Eq. 

In the light of these considerations, we can more easily 
interpret the results of the DFA for the two-dimensional 
Hubbard model, which are presented in Figs. [SJ [HI [IHl 
Specifically, we start the analysis of the two dimensional 
case, by evaluating the effects of the Moriya A correction 
for the DFA results computed for the half-filled Hubbard 
model with ?7 = 4i at a temperature (/3 — 17) slightly 
above the corresponding T/v in DMFT. 

In the first/third row of Fig. [51 we show the DFA self- 
energy and spectral function at the Fermi surface (FS) 
at the nodal [q = (^, |)]/antinodal [q — (tt, 0)] points 
computed without Moriya correction. One can clearly 
observe that, in contrast to the three dimensional case, 
the DFA spectra qualitatively differ from the original 
DMFT one because (i) a pseudogap appears at low fre- 



7 



Res. 



Im 






U=D 


DMFT -- 


■q=(it;2,7i/2) 


J, Qpy^ — 







Re I|, 



Im 



DMFT . 

X-DrA - 



U=D 

q=(n,0) 



aL-a 



2 2 
co/D 



-0.5 




(o/D 



0.5 



FIG. 8: (Color online) DFA results for the half-filled two- 
dimensional Hubbard model a,t U = D — 4t, (3 = 17 (just 
slightly above T^^^'^) computed without (first and third 
row; black dotted line) and with (second and fourth row; 
blue solid line) the Moriya A correction, and compared with 
DMFT (grey dashed line) . The DFA calculations in the non- 
selfconsistent scheme show a clear pseudogap opening for the 
k points of the non-interacting FS, but more pronounced in 
the antinodal direction. Within the A-corrected scheme, one 
can still notice a pseudogap opening but only in the antinodal 
direction, while in the nodal direction a strongly damped QP 
appears. 



quencies and (ii) the spectra are markedly anisotropic in 
the nodal/antinodal direction, as the observed pseudogap 
is evidently more pronounced at the antinodal points. 

As discussed above, the Moriya A correction is how- 
ever expected to be much more important in the two- 
dimensional than in the three dimensional case . This 
is confirmed by the results shown in the second/fourth 
row of Fig. [H In these panels the reduction of the non- 
local effects due to the inclusion of the Moriya A correc- 
tion in DFA is evident. It is important noticing, how- 
ever, that although the distance in the phase-diagram 
from the actual antiferromagnetic transition (occurring 
at T = for d = 2) is considerably larger than for 
d — 3, non-local correlation effects are nonetheless ex- 
tremely strong. Turning to the details, we still observe a 
remarkable anisotropy in the DFA spectra after the inclu- 
sion of Moriya correction, with a strongly renormalized 
QP peak in the nodal direction and a rather clear pseudo- 
gap in the antinodal direction. The results of DFA (im- 
plemented with the Moriya correction) indicate, hence, 
that for the two-dimensional system at half-filling anti- 
ferromagnetic fluctuation effects predominate in a wide 
region of the phase diagram, determining the onset of an 
anisotropic pseudogap in the spectra also for considerably 
high temperatures, qualitatively similar to that observed 
in underdoped cuprates.— 

The inclusion of the Moriya correction in DFA allows 
us to extend our analysis to the low-temperature regime 
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FIG. 9: (Color online) Temperature evolution of the DFA 
results for the half-filled two-dimensional Hubbard model at 
U = D = 4Hn the nodal direction, computed with the Moriya 
A correction, and compared with DMFT. A clear pseudogap 
emerges at the lowest temperature (/3 = 60), similarly to the 
results of the non-self consistent scheme in the proximity of 
T^^^^ (see Fig.ia. 



T < T^^^'^. In particular, we are interested to study 
the evolution of the spectral function when the temper- 
ature is considerably reduced compared to the DMFT 
value T^^^^. In Figs. [Hand [TO] we report the DFA caF 
culation for the self-energy and the spectral function for 
the same case considered above {U = At, half-filling) for 
three different decreasing temperatures (/? — 17, shown 
already in Fig. [SJ /? = 25 and (3 — 60) in the nodal and 
antinodal direction, respectively. First, we note that the 
anisotropy in the self-energy and the spectra remains vis- 
ible at all temperatures. In addition, a marked tendency 
towards a completely gapped spectrum can be seen at 
the lowest temperature: At lowest temperature (/3 = 60) 
a pseudogap appears also in the nodal direction, while 
the pseudogap already present in the antinodal direction 
becomes remarkably more profound. At this tempera- 
ture, therefore, the anisotropy is reduced in comparison 
to the higher T cases- due to the strong depletion of spec- 
tral weight at a; = 0. This results can be understood in 
terms of the closer proximity to the antiferromagnetic 
instability at T — 0, and is consistent with the marked 
pseudogap visible in the k-integrated spectral function 
obtained by means of cluster DMFT in Ref. [TtI . 

It is worth noticing, however, that the temperature 
evolution towards the formation of a fully gapped spec- 
trum at T — + does not appear to be completely 
monotonous. The effects of the non-local fluctuations 
seems to be slightly weaker in the DFA results for /3 — 25 
(second row in Figs. [Ql fTO)) . than for — 17 (first 
row). More specifically, this is visible in the slightly more 
Fermi-liquid-like behavior of the real and imaginary part 
of the self-energies at /3 = 25 in comparison to (3 = 17. 

A possible interpretation of this specific feature of 
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FIG. 10: (Color online) Same as in Fig. [9] in the antinodal 
direction. As expected, a very pronounced pseudogap char- 
acterizes the lowest temperature results. The behavior of the 
spectral functions is not completely monotonous, as the pseu- 
dogap seems to disappear at /3 = 25. At all temperatures, 
however, the pseudogap features are always more marked in 
the antinodal than in the nodal direction. 
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FIG. 11: Color online) DFA results with A corrections for the 
half- filled two-dimensional Hubbard model at ?7 = 2D = 8t, 
l3 = 40 compared with the corresponding DMFT ones. 



our results is to relate the non-monotonous tempera- 
ture evolution in the DFA spectral function to a com- 
petition between non-local and local mechanisms capa- 
ble of destroying coherent excitations: (i) The (non- 
local) antiferromagnetic fluctuations, which become less 
pronounced with increasing T, making the system more 
metallic (xq,o = 8.9 • 10^, 39.26, and 13.28, for (3 = 60, 
25, and 17, respectively); and at the same time (ii) the 
thermal loss of coherence, which is at the origin of so- 
called crossover region in the (purely local) DMFT and 
reflects increasing values of the quasiparticle damping 
(7 = -lmE(O) = 0.009,0.021,0.034, for the three con- 
sidered temperatures respectively) . The relevance of the 
interplay between these two mechanisms is an interesting 



issue raised by our DFA results. It might also be related 
to a similar non-monotonous trend in the cluster DMFT 
phase diagram reported by Park et alM. 

The DFA results at stronger interaction {U = 2D and 
f3 — 40) are presented in Fig. [TlJ At the considered 
low temperature the local DMFT spectral functions have 
peaky structure, because we solve the impurity problem 
of DMFT by means of exact diagonalization (ED), which 
treats only finite number of sites. Note, however, the 
DFA spectral functions are continuous due to momenta- 
and frequency sums in the Eq. ([8]), even though ED 
is employed as an impurity solver. The nonlocal spectral 
functions show the splitting of the quasiparticle peak due 
to magnetic correlations, which is similar to the structure 
(d) in Fig. [3] discussed in Sect. IV^^. Closer to the Mott 
transition (i.e. at even stronger U) we also expect the 
formation of the structures (a)-(c) of Fig. [3] 

The presented results demonstrate that the DFA ap- 
proach -with the inclusion of the Moriya corrections- al- 
lows for a non trivial analysis of the effects of long-range 
spatial correlations in every region of the phase diagrams 
of strongly interacting fermionic systems both in two and 
three dimensions. 



VIII. k-RESOLVED SPECTRAL FUNCTIONS IN 
TWO DIMENSIONS 

Let us now calculate the k-dependence of the spectral 
functions in the directions of high symmetry, as can be 
observed in angular resolved photocmission spectroscopy 
(ARPES). It is worthwhile remarking that, in contrast 
to the cluster extensions of DMFT, this does not require 
any kind of interpolation in k-space: Due to the diagram- 
matic nature of the DFA, the spectra for every chosen k 
point in the first Brillouin zone are easily computed via 
Eq. ^. 

Here, in Fig. [121 we present DFA results with Moriya 
A correction for the same case previously considered in 
Fig. [8] (second and fourth row). As it is often done, we 
consider two different k-paths along the Brillouin zone, 
the first one along the nodal direction [(0,0) (7r,7r), 
left panel] and the second one right at the border of the 
Brillouin zone, crossing the antinodal point at the FS 
[(7r,7r) — > (tt, 0), right panel]. 

Our analysis of the k-resolved DFA results allows us to 
appreciate the evolution of the main features of the DFA 
spectral functions. Specifically, we observe that for the 
points most far away from the FS, the spectral functions 
display similar features in the two cases: A relatively nar- 
row peak separated from a broader maximum at higher 
energies, which represents the incoherent processes build- 
ing up the (upper) Hubbard band. When proceeding in 
the direction of the FS, as expected, the narrow peak 
moves towards the Fermi energy, while the broad feature 
becomes less pronounced. A qualitative difference be- 
tween the two selected paths emerges only in the vicinity 
of the FS: The shift of the narrow peak down to zero en- 
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FIG. 12: (Color online) k-resolved DFA spectra along the 
nodal (left) and antinodal (right) direction for the half-filled 
two-dimensional Hubbard model at f/ = D = it and /3 = 17, 
calculated with the Moriya A correction. 



ergy is frozen along the second path, consistent with the 
opening of the anisotropic pseudogap in the antinodal 
direction, while it continues to shift down to the Fermi 
level in the nodal directions. 

It is also worth noticing the occurrence in both cases of 
a slight broadening of the narrow peak while approach- 
ing the FS. This trend, which is markedly different from 
any FL expectation, could be understood in terms of the 
maximum of Im S(k, tu) appearing at zero frequency (see 
again Fig. [8|) for both directions. The enhanced value 
of ImS(k, w) at low frequencies, which is ultimately re- 
sponsible for the opening of the pseudogap starting at 
the antinodal points, determines a loss of coherence and, 
hence, the observed broadening of the peak, while it 
moves closer to the FS. 



IX. CONCLUSION 

Based on the representation of the nonlocal self energy 
which considers the effect of the bare Coulomb interac- 
tion and charge (spin) fluctuations, we have extended 
the recently introduced dynamical vertex approximation 



(DFA) by including a Moriyaesque A correction to the 
local vertex in Section |Vl The value of A is determined 
from the sum rule which relates w-integrated self en- 
ergy and occupation and allows for a proper reduction 
of the DMFT Ncel temperature, in two dimensions even 
to Tn = so that the Mcrmin- Wagner theorem is ful- 
filled. This correction is therefore particularly important 
for two dimensions, where spin fluctuations are especially 
strong. Without the Moriya A correction, a much more 
involved self-consistent solution of the DFA equations 
would be necessary to yield similar results. 

The method we have introduced here allows for a treat- 
ment of non local long-range spatial correlation in finite 
dimensional systems. In three dimensions, pronounced 
effects of non-local spin fluctuations are found only close 
to the antiferromagnctic phase transition. This is in 
contrast to the two dimensional case where antiferro- 
magnetic fluctuations completely reshuffle the spectrum, 
also far away from the antiferromagnctic phase transi- 
tion at Tn — 0, leading eventually to the formation of a 
pseudogap. Qualitatively, the spectral functions can be 
understood by means of the analytical formula for the 
self energy proposed in Section |TVl Calculating several 
DFA self energies along the high symmetry lines of the 
Brillouin zone, we obtain the momentum dependence of 
the spectral functions, which could be directly compared 
with the ARPES data. 

DFA can serve as a very promising method for fu- 
ture studies of the Hubbard model at non-integer fill- 
ing, in particular in the vicinity of the antiferromagnctic 
quantum critical point i^i^^ A further important devel- 
opment would be also the generalization of the method 
to the multi-orbital case, to analyze the effects of non lo- 
cal correlations beyond DMFT in realistic bandstructure 
calculationsi^ 
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